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^_ ■ Abstract 

de Sitter space is known to have a thermal character. This can be best seen by an 
Unruh-DeWitt detector which stays in the Poincare patch and interacts with a scalar 

IT) ! 

field in the Bunch-Davies vacuum. However, since the Bunch-Davies vacuum is the 

ground state only at the infinite past, if the scalar field starts in the ground state at 

a finite past, an Unruh-DeWitt detector then will feel as if it is in a medium that is 

^ \ not in thermodynamic equilibrium and that undergoes a relaxation to the equilibrium 

5—i < 

corresponding to the Bunch-Davies vacuum. In this paper, we first develop a general 
framework to treat such circumstances and write down the master equation which 
completely describes the finite time evolution of the density matrix of an Unruh-DeWitt 
detector in arbitrary background geometry. We then apply this framework to an ideal 
detector in de Sitter space which can get adjusted to its environment instantaneously, 
and show that the density distribution of the detector certainly exhibits a relaxation to 
the Gibbs distribution with the universal relaxation time of half the curvature radius 
of de Sitter space. This relaxation time gives the first example of such quantities that 
are related to nonequilibrium dynamics intrinsic to de Sitter space. 
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1. Introduction 



The concept of particle is known to depend on observers. Even in the Poincare-invariant 
Minkowski vacuum, an observer with constant acceleration sees a thermal particle spectrum 
PP . The Unruh-DeWitt detector [TJ 12] is invented as a tool of thought experiment to give an 
intuitive understanding of such thermal character of spacetime (see also [31 H] and references 
therein). This is a detector weakly interacting with a matter quantum field in a certain 
vacuum state, and one can study the thermal character through the density distribution of 
the detector. 

Various spacetimes have been examined with the Unruh-DeWitt detector, including de 
Sitter space. There are a family of de Sitter-invariant vacua (called the a-vacua) for a free 
scalar field in de Sitter space 0E]- Among them, the most natural vacuum, satisfying the 
Hadamard condition, is the Bunch-Davies vacuum (or the Euclidean vacuum) [7], which 
actually exhibits a thermal property. In fact, an Unruh-DeWitt detector staying at x = 
in the Poincare patch 

ds 2 = f ~ dr l 2 + dx2 ( oo<r]<0 ) (i.i) 

responds as if it is in a thermal bath with temperature T = l/2ir£ , where I is the curvature 
radius of de Sitter space [3j EJ E] (see also [10J ITT]). 

In our previous paper [12] , it is shown that the Bunch-Davies vacuum corresponds to the 
configuration where the scalar field starts in the ground state at the infinite past (77 — > —00). 
If instead the scalar field starts in the ground state at a finite past (say, at 770) , then the de 
Sitter invariance should be broken explicitly, and we expect that an Unruh-DeWitt detector 
feels that the surrounding medium is not in thermodynamic equilibrium and undergoes a 
relaxation to the equilibrium corresponding to the Bunch-Davies vacuum. 

The main purpose of this paper is to develop a machinery for describing such nonequilib- 
rium dynamics and to calculate quantities directly related to the nonequilibrium properties 
intrinsic to de Sitter space. For this, we first develop a general framework to describe an 
Unruh-DeWitt detector in arbitrary background geometry and derive the master equation 
which completely determines the finite time evolution of the density matrix of the detector. 
We then apply this framework to an ideal detector in de Sitter space which can get adjusted 
to its environment instantaneously. We show that the density distribution of such detector 
exhibits a relaxation to the Gibbs distribution with the universal relaxation time £/2 mea- 
sured in the proper time of the detector. This relaxation time gives the first example of such 
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quantities that are related to nonequilibrium dynamics intrinsic to de Sitter space. 

This paper is organized as follows. In section [21 adopting the method of the projection 
operator [13], we first derive the master equation which describes the time evolution of the 
density matrix of the detector. Then, after justifying a Markovian approximation, we derive 
a simplified form of the master equation which enables us to analytically study the relaxation 
behavior of the density distribution. In section [3j we perform an analytical calculation of 
the transition rate matrix, expanding it around rjo = — oo . In section HJ we consider an ideal 
detector with two energy levels, and show that it is always thermalized with the relaxation 
time £/2. Section \5\ is devoted to discussions and conclusion. We collect miscellaneous 
formulas in appendices. 

2. Master equation for the density matrix of an Unruh- 
DeWitt detector 



field 4>(x) of mass m . We assume t 
it moves along a classical trajector^ 



2.1. Setup 

We consider an Unruh-DeWitt detector in d-dimensional spacetime with background metric 
ds 2 = 9[iv{x) dx M dx u (fi, v = 0, 1, . . . , d — 1), the detector being interacting with a scalar 

Jiat the detector has a sufficiently large mass so that 
x^{r) = (t(r), a;(r)) , where r is a proper time of the 
trajectory (see figured]). 

For a quantum mechanical description of the system, we introduce the Hilbert space T-L tot 
which is the tensor product of those of the detector and the field, 

U tot = H d . (2.1) 

The total Hamiltonian then takes the following form in the Schrodinger picture: 

H tot (t) = H d ® 1 + 1 ® H+(t) + V(t) . (2.2) 

Here, H d is the Hamiltonian of the detector associated with its proper time, and we assume 
that H d does not depend on r (or on t), denoting its eigenstates and eigenvalues by \m) 
and E m , respectively: H d \m) = E m \m) . H^(t) is the Hamiltonian of the free scalar field 

1 We also write the trajectory as x^(t) = (t 7 x(t)) using the functional relation t = t(r) or r = r(£) with 
dt/dr > 0. 
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Figure 1: The trajectory of an Unruh-DeWitt detector, which starts interacting 
with a scalar field 4>(x) at time t\ . The scalar field will be assumed to be in the 
ground state at an early time to • 



and may depend on time explicitly through an explicit time dependence of the metric. V(t) 
stands for the interaction between the detector and the field, and we assume that it suddenly 
starts at time t\ in the form of monopole interactional 

V(t) = A n ^ ® <f>(x{t)) 9(t - h) . (2.3) 

Here, A is a dimensionless coupling constant, x(t) represents the position of the detector at 
time t (see footnote [T]) , and \i is an operator acting on "H d , which we again assume to be 
time independent. The time evolution of the total density matrix p tot (t) is given by 

p tot ^ = jjtot ^ ^ p tot ^ ptot ( t ,t')]~\ (2.4) 

where U tot (t,t') is the time evolution operator in the Schrodinger picture, 

U tot (t,t') = Texp(-i^dt'# tot (t')) . (2.5) 

We consider a situation where one can only measure observables associated with the 
detector, such as the matrix elements \i mn = (m\n\n). Then the maximum information one 
can get from the system is the reduced density matrix p(t) which is defined as the partial 
trace of the total density matrix p tot (t) over Jf| 

p(t)EETr p tot (t). (2.6) 



2 Recall that this is in the Schrodinger picture. We later shall take an average over the start-up time t\ 
since one usually needs a certain period of time to specify the initial distribution of the system. 

3 More precisely, Tr^ is a linear map from an operator of the form O = ^2 O d <£> Of Y , which acts on H tot , 

a 

to the operator Tr O = J2 0« x £ <-0|O^|V> , which acts on Ti d , where £ \ip){ip\ = 1 . 



Since the dynamics of p tot (t) is completely fixed, the time evolution of p(t) should be uniquely 
determined once one specifies the form of interaction and the initial condition for p tot (t) . 
Since we assume that there had been no interaction between the detector and the field before 
time t\ , we may let the total density matrix take the following factorized form at t\ : 

p^fa) = p d (ti) ® p^i) = p{h) <g> /(tO . (2.7) 

The latter equality can be easily seen by noting p{t\) = Tr^, p tot (tO = P d (^0 ■ 



2.2. Master equation 

The time evolution of p(t) can be best analyzed if we go over to the interaction picture by 
decomposing the total Hamiltonian to 

H tot (t) = H^it) + V{t) , (2.8) 

m ot (t) = H d ® 1 + 1 <g> H+(t) (2.9) 

dt 

and treat V(t) as a perturbation. The time evolution operator U tot (t, ti) is then decomposed 
to the unperturbed and perturbed parts as 

f/ tot (M0 = c/ tot (M0 Uf\tM) (2.10) 

with 

C/ O tot (t,ti) = Texp(-i jTat'^Ct')) = e - iHi < T ^®Te- l & d * H *V ) , (2.11) 

C/) ot (t,t0 = Texp(-i jf^'V^f')) • (2.12) 



Here, V/(t) is defined by 



V^(*) = [^(Mi)]~V(i) ^(Mi) 



= A ^ /i 7 (r) <8> 0j(z(r)) 0(t - t0 , (2.13) 

where pi(r) is given by 

pj(T) = e^-f-O^-ifMr-n) = ^ e i { E m -E n){T - T1 ) ^ | m)(n | ? (2 w) 

and the operator 

h{x(r)) = = [Te-^ d ^ ( '' ) ]-V(^W)Te" i ^ d ^ ( '' ) (2.15) 
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satisfies the free Klein-Gordon equation associated with the metric ds 2 = g^ix) dx M &x u . 
Accordingly, the density matrix in the interaction picture is given by 

P f\t) = [uF&tiT^WuF&ti) = ^ ot (*,*i)p tot (*i) [fTCMi)]" 1 , ( 2 - 16 ) 

and satisfies the von Neumann equation 

±pf\t) = - i [Vj(t),pr(t)] = -i ad y/(i) pf\t) . (2.17) 

The above von Neumann equation can be rewritten to an equation involving only pi(t) = 
Tr^pj ot (i) by adopting a well known procedure called the projection operator method (see, 
e.g., PS])- We first introduce a linear operator V: End"H tot — > End"H tot which acts linearly 
on elements belonging to End"H tot (the set of linear operators acting on "H tot ) and has the 
form 

V : O ^ VO = (Tr O) (O e End H tot ) . (2. 18) 

Here, e End IH/ can be any operator acting on 'W' as long as it satisfies 

Ti>X* = l (2.19) 

and does not depend on time. From (12.191) one can easily see that V is a projection operator, 
V 2 = V. For O = pf{t), we obtain 

Vpf\t) = (Tr^pf(t)) ® X* = Pl {t) ® X* . (2.20) 

We further introduce Q = 1 — V , which is also a projection operator, Q 2 = Q , and satisfies 
VQ = = QV. 

Since P/ ot (t) satisfies the time evolution equation (I2.17p . Vp^it) and Qp^it) satisfy the 
differential equations (note that V + Q = 1) 

j t Vpf\t) = -iPui Vl ® V P T(t) - iVad Vl{t) Qpf\t) , (2.21) 

j t Q P T(t) = -iQad Vl(t) V P T(t) - iQad Vl(t) Qpf\t) . (2.22) 

The solution of (I2.22p is given by 

Qp\ ot (t) = -i / dt' Te' 1 $ dt " Sad ^ (t ") Q ad Vj(t0 Ppf^') + Te _i d '' Qad W> Qpf 1 ^) , 
Jti 

(2.23) 
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and substituting this to (I2T211 and using (&/ &t)V p\ ot {t) = (dpj(t)/di) <8> , we obtain 

-^ad V/W / dt'Te- i ti dt " QaA W')Qeid Vlit , ) Vp¥ t (t'). (2.24) 

This equation can be simplified by noticing that can be chosen arbitrarily without 
changing the time evolution of pi(t) . We here set 

X* = pfih) , (2.25) 

which certainly satisfies the condition (12.191) . Then, due to the initial condition (12. 7p . we 
have 

VpTih) = pT{h) , QpTih) = , (2.26) 

and thus the second term on the right-hand side of (12.241) vanishes. Furthermore, since the 
interaction has the factorized form (I2.13p . we obtain 

T>(ad y/(t) Vpfit)) = A ^ D*i(t), Pl (t)] Tr (0 7 (x(r))pf (t x )) , (2.27) 
which vanishes when (as we assume hereafter) the one-point function of the scalar field is 



Ti^ I (x(T))pt(t 1 ))=0. (2.28) 

We thus find that the first term in (12.241) also vanishes, 

V ad Vj(t) VpT(t) = Tr (adv^) VpT(t)) ® X* = . (2.29) 

The last equation f!2T2"9l shows that Qad Vl{t) = ad Vl{t) Vpf\t) = ad V/(t) (p 7 (i) ® 

Pj{t\j) , and thus, (I2.24p becomes 



dpj(t) 
dt 



pf(ti) = -Vad m) f dt'Te-'ti* 1 "^^ ad Vj(t0 (p,(f) ® pf(*i)) . (2.30) 



This is the master equation which with the initial condition (12 .7p and the assumption (I2.28P 
completely determines the time evolution of the reduced density matrix pj{t) . 
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2.3. Approximation of the master equation 

We expand the right-hand side of (12.301) to the second order of the perturbation theory to 
obtain 



dpi(t) 
dt 



Tr^ (ad Vl (t) £ dt'ad Vl{tl) (p/(0 ® + ^) • ( 2 -31) 



This can be further rewritten in terms of the proper time r to the following form (denoting 
p(t(r)) by p(r) and derivatives with respect to r by dots): 

Pi(r) 

= - A 2 rdr'Tr4^(r)®0/(x(r)), [p^r') ® 0,(x(r')), P i{t') ® pj(n)]] + 0(A 3 ) 

= A 2 f dr'(W), wMMO] G+(x(r'),*(r)) - ^(r), M^VCr')] (x(t), x{t')) 
+ 0(A 3 ) , (2.32) 
where 

G\{x,x') = Tr^Mx) M^)Pi(n)) (2.33) 

is the Wightman function of the free scalar field with the density matrix = pf(ri) . 

Since the reduced density matrix in the interaction picture, pi(t), is related to that in 
the Schrodinger picture, pit), as 

Pl ( T ) = e iRi - (r ~ Tl) P (t) e~ iRi - (r ~ Tl) , (2.34) 

we can rewrite (I2.32p to the following form in the Schrodinger picture: 

p(r) + i [H\ p(r)] = A 2 f dr> [p p(r >) ^ G + ^ x(r)) 

+ e~ iHd - W p p(r') p e iHd - ^ G + x (x(r), x(r')) 

_ e -\H A - (t-t') n p iH d - (r-r') 



p(r')pe iHa -^pG+(x(T , ),x(r)) 



- p e-' lHd - p p(r') e iHd - ^ G x (x(r), x(r'))] 
+ C(A 3 ) , (2.35) 
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which has the following expression in terms of matrix elements: 



p mn {r) + i(E m - E 

n) Pmn \T} 

k,i Jti 

+ e-HZn-mr-*)^ ^ pkl{j >) G + (x(r), x(r')) 

- e~ l ^- E ^'-^ m urn Pmk(r') G + x (x(r'), x{r)) 

- e-^-^^'W/^ Pm(r') G x (x(t),x(t')) } + 0(\ 3 ) . (2.36) 

The integro-differential equation (12.361) can be further simplified as follows. Since the 
Wightman function G x [x(r), ^(t')) in the integral is singular at r' = r and decreases 
exponentially for the large separation of r and r' the main contributions to the integral 
should come only from the region r' ~ r . This implies that the memory effect in the equation 
is highly suppressed, and thus we may replace Pfcz( r ') m the integral by its boundary values 
Pkiij) to a good accuracy, assuming that p(r') slowly changes. We then have 

Pmnij) + \{E m - E n )p mn {r) 

- f ^'[^ {En ~ Ek){T '~ T) PmkPinPki{T)G x (x(r'),x(r)) 

k.l 



+ e -i( £m - £l )( T -r') tafem(T) G+(x(t),x(t')) 
Atfei Pin Pmfc(r) G^(x(r'),x(r)) 



_ e -i( Si -S m )( r '-r) 
_ -i(E k -E n )(r-T') 



Pmk Pm Pin(r) G x (x(t),x(t')) ] . (2.37) 



If we 



"urther assume that the off-diagonal elements of p(r) almost vanish in the region 
t' ~ r.f| then (I2.37P becomes 



^ ^2 \^rnk\ 2 / dr' 
k^m 

x [e-^^-^^VttW G + x (x(r'), x(r)) + e -i(M*)(r-r ') pfefc(r) G + ( x(r); x(r ')) 



4 We see in Appendix [B]that for a scalar field in de Sitter space the Wightman function certainly exhibits 
this property if the mass is large enough. In general, there can be a case where the Wightman function has 
a long tail and one needs to take account of memory effects carefully. We do not deal with such cases in the 
present paper. 

5 This assumption may seem to be too strong, but we will see in section 2] that we can always set the 
off-diagonal elements to zero if the detector is a two-level system and \i is off-diagonal. 
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- e-^- E ^'-r) pmm{r) Q+x ( x (t>), x {t)) - e -i(*fc-*»)(r-T ') pmm (r) G + x (x(r), x(r')) ] 
= ^2 [ w mk{r, Ti)p kk {T) - w km (r, ri)p mm (r)] . (2.38) 

Here, we have introduced the transition rate matrix 

Wmfc(7,Ti) = A 2 \fi m k\ 2 F{E m - E k ;r, ri) , (2.39) 

where 

F{AE\ r, n) = / T dr' [ e - iAB(T '- T) G+ (x(r'), z(r)) + e - iAS (^') G+ (x(r), x(r'))] . (2.40) 

Equation ( 12.381) has the standard form of the master equation. If, in particular, J req (A£ l ) e 
lim ^(AE; r; ri) satisfies the relation^ 

r— ri— S-oo 

^(AE) ^ 
J C( i(-A£) 

the transition rate matrix satisfies the relation 

W m k(r; Tj) r-7i->oo e -^( Em -E k )_ (2.42) 
Wfcm(T";Ti) 

Then, when the transition rate matrix w mk (T\T\) is irreducible (i.e., if /i does not take a 
block diagonal form with respect to the basis \m)), the distribution of the detector in equi- 
librium, p eq , can be determined by the detailed balance condition lim w mk (r\Ti) p e ^ k = 
lim w km (r; n) p„ m , and we obtain 

r — 71— >oo 



p — pEm ^ 

n 

which is nothing but the Gibbs distribution at temperature 1//3. 

We close this section by making a comment on the relationship between our formalism 
and the literature. One can easily show that the transition rate matrix w mk , (12.391) . is the 
derivative of the following expression with respect to the proper time r: 

X 2 \fi mk \ 2 F{E m -E k ]T,T 1 ), (2.44) 
F(AE; r, n) = f dr' C fa» e -^V -*") £+ (x(t'),x(t")) . (2.45) 



6 We will see that this condition is indeed satisfied with (3 = 2ir£ when p^{t\) corresponds to the Bunch- 
Davies vacuum in dc Sitter space. 
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In the literature (e.g. [5]), one often considers a process starting from the initial state |/c)®|a) 
at time t\ (usually taken to be the infinite past) to the final state \m) £g> \/3) at time t, 
and sums over the final states of the field, \/3) . The transition probability has the same 
form as (I2.44p if we set = p^(ti) = | q;) (ck | for which the Wightman function becomes 
G|-(i, x') = (a\ 4>i(x) <f>i(x') \a) . We thus again see that (12.39}) represents the transition 
probability per unit proper time of the detector. J 7 in (I2.45P is often called the response 
function (see, e.g., [3]). As we have seen, J 7 or its derivative J- does not depend on details of 
the detector and can be thoroughly determined by the Wightman function of the free scalar 
field. 



3. Unruh-DeWitt detector in de Sitter space 

In this section we consider an Unruh-DeWitt detector in d- dimensional de Sitter space, which 
is weakly interacting with a massive scalar field 4>(x) of mass m . We exclusively consider 
the Poincare patch denoting the time variable by r\ , 

ds 2 = e 2 ~ dr ) 2 + dx2 (_ 00<r? <o), (3.i) 

and set the classical trajectory of the detector to be the geodesic 



.r" 



(r) = (-£e-^,0). (3.2) 



We denote by 771 the time when the detector starts the interaction with the field <fi{x) = 
(j)(r],x), which has the following form in the interaction picture (see (12.131) ): 

Vrfo) = A N (r) g, 0/ ( x(r) ) g {rj _ m) . (3.3) 

We assume that the scalar field is in the instantaneous ground state |O w ,77o) a t an early 
time 7/0 (see figure. [2]). We introduce the state \0 VO ) = 10^,771) as an interaction picture 
state (or as a Heisenberg picture state of free field theory) which is obtained by applying the 
free-field time evolution operator from r] to r]i : (0^) = Texp(— i f v *dr]' if ^(7/)) |0 % , 77 ) E 

We now set = lO^XO^J . Note that the one-point function certainly vanishes, 
Tr0 0(x)X^ = (0^ I 4>{x) \Q m ) = 0. Then, with the approximations made in subsection 



7 See [12] for a detailed discussion on defining such time-dependent ground states for a free scalar field 
in a curved spacctimc. 
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Or, > = lO^o, r?i) 

Figure 2: The trajectory of an Unruh-DeWitt detector in the Penrose diagram 
of de Sitter space. The shaded region corresponds to the Poincare patch, and 
the dashed line represents the future event horizon for the detector. The scalar 
field is in the ground state at time r] . 

\2.3\ we have the master equation of the form (see (I2.38p - fl2.40p ) 

) = ^ [wmk{r, Ti, r]o)pkk(r) - w km (r, n; r/ ) p mm (r)] (3.4) 

with the transition rate matrix 

w mk (r, n; rjo) = X 2 \pmk\ 2 F{E m - E k ; r, n; T] ) . (3.5) 
Here, the derivative of the response function takes the form 

F(AE;T,n; Vo ) = f dr'[e^ E ^'^ G + (x(r'), x(r); + e ~ iAE ^ G + (x(r), x(r'); %)] 
J n 

(3.6) 

with the Wightman function 

G + (x,x'; m ) = (OJM^M^K) ■ ( 3 - 7 ) 

In this section, we first give a brief review on how to obtain G + (x, x'; t]q) for finite rjo closely 
following [T2] and then calculate ^(AE; r, T\, r/ )E 

3.1. Wightman function 

In the Poincare patch, there is a translational invariance in the spatial directions, and we 
assume that all spatial directions are compactified with radius L/2ir (which will be taken 
to infinity in the end). Then the wave vector k takes discrete values k = [2i\jV)n with 

8 In the following, we set the curvature radius I = 1. 
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n G Z d 1 . We denote k > (or k < 0) if the first nonvanishing element of the vector 
fc = k 2 , ■ ■ ■ , k d _i) is positive (or negative). Now, we consider the mode expansion 

<f)(x) = (f>(r), x) = ^2 faM Y k,a{x) , (3.8) 

fc>0 a 

where the mode functions \Y k>a (x)} are defined by 

k = 0: Y k=0 , a=l = ^= (V = L d - 1 ), (3.9) 
v V 



k>0: Y k>a= i(x) = y ^ cos(fc • x) , ife ;0=2 (cc) = y y sm(k ■ x) . (3.10) 

Here, {Y kja (x)} is a complete set of (real-valued) eigenfunctions of the spatial Laplacian 
Arf-i = 5^i=i ^ 2 and satisfies the orthonormal relations 

[d d - 1 xY kta (x)Y kl>a ,(x) = 5 k>kl 5 aja ,. (3.11) 

With the mode expansion ( 13. 81) . the Hamiltonian (in the Heisenberg picture of free field 
theory) is expressed as the sum of the Hamiltonians for independent harmonic oscillators: 

H (v) = Wk(v) a k,a(v) + const - > ( 3 - 12 ) 



fc,a 



<*kM = \j P ^^ k ^ (f>kM + i\j 2 p (rj)u k {rj) nk ' a ^ ' n kM = p(v) <f>k,a(v) , ( 3 - 13 ) 



= (_^-(d-2) , Wfc ( 77 ) = <N / m 2(_^)-2 + A .2 (fc=|fc|). (3.14) 

The instantaneous ground state (0^) is then expressed as a k}a (Vo) \0 Vo ) = . 
The free Klein-Gordon equation takes the fornix 

rf 4> k ,a(v) ~ (d ~ 2) j> k ,a{v) + 0V + ™ 2 ) ^(77) = , (3.15) 

and we denote a pair of independent c- number solutions by {f(rj), g(r))} . One can easily 
check that their weighted Wronskian W p [f, g] = p{rf) [f(rj) g(rj) — f{rf) g{rj)\ does not depend 
on r\ . By expanding the operator <j) kt a{v) using the annihilation and creation operators 
a k , a (Vo) and a^ a (%) as 

<t>k,a(v) = Vkiv, Vo) akAvo) + <Pk(v, Vo) a lJvo) , (3.16) 



9 Note that, in calculating the in-in propagator in the Poincare patch, we do not need to add — ie to m 2 
because there is no subtlety in specifying the ground state wave function for this case, as can be seen from 
calculations performed in subsection 4.2.1 of [12]. 
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the ground state wave function tpkiv'iVo) a t time 77 is obtained in the following form 

1 



[vofiv) -uog{v)] 



(3.17) 



u = p(r] ) [/(770) + iuJkivo) f(Vo) ], v = p(rj ) [g(rj ) + iu k (r] ) g(rj ) ] . (3.18) 
The wave functions can be shown to be independent of the choice of the pair {f(rj), g(fj)} 

im 



We now choose clS db set of independent solutions 

f(v) = (-V)^ U-kv) , 9(V) = (-V)^ N u (-k V ) , (3.19) 
where J u (x) and N u (x) are the Bessel and Neumann functions, respectively, and 

d- 1\ 

(3.20) 



v = < 



d- 1\ 2 



< m < 



1 a / m z 



d-l\2 / d-1 

m > 



2 / V - 2 

For this pair of solutions, the weighted Wronskian is given by W p [f, g] = —2/n, and, uo and 
Vo in (13.181) take the forms 

■d- 1 



u 

^0 



-w 2 

-Vo) 2 



2 

d- 1 



+ u+ iuj k (r]o)Vo) Ju(-krjo) + k rj J 1+u (-k rj ) 



+ z/+ ia; fc (?7o)77o) N v (-kr) ) + kr] N 1+u (-krj ) 



(3.21) 
(3.22) 



The wave function is then given by 



7T 



■v) 



d-1 
2 



, d-2 



2i/2o; fe (77o, 

from which the Wightman function for each mode is given by 



[vo Ju(-k 77) - u N u (-k 77)] 



(3.23) 



Gt{r],r]';r) ) 



(0,0 


<l>k,a(v) 


0fc,a(V) 


W ) 


(0 W 


0,o) 



(3.24) 



Summing this over k , we obtain the Wightman function in spacetime as 
G+(x,x'; Vo ) = ^ Gt(v,v';Vo) Yl Y U X ) Y kA x ') 



k>0 



d^k 



cos [k ■ (x - x')] ip k (r]; rj ) ip* k (rj'; rj ) 



(3.25) 
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where we have taken the limit L — >■ oo in the last equality. Integration over angular variables 
can be easily performed [12], and we have the expression 



G + (x,x';r)o) 



1 



(2ir) 2 \x — x | 



2 JO 



dkk"^ Jd^i (k\x- x'\) (p k (rj; rjo) fl(f]'l Vo) 



(3.26) 



It is shown in [12] that the Wightman function becomes that for the Bunch-Davies vacuum 
in the limit r] — >• — oo : 



lim G + (x, x'; rjo) 



4 



(2tt) 2 \x — x'\ 



2 Jo 



dkk^ J^z(k\x- x'\) H^i-kr]) Hl 2 \-kr]') . 



r (¥ + ^) r (¥~^) p (d-l . r , d-1 d.l-ti 
(4vr) rf / 2 r(f) 2 l( v 2 + ' 2 '2' 2 

— ^bd( 3; ' x ) ) 



(3.27) 



where u = — [(77 — 77' — iO) 2 — (a? — a;') 2 ] /(2 77 rj') and 2-^1 ( a , c; z) is the hypergeometric 
function. In the next subsection, we expand the wave function (fkiv'iVo) around rjo = —00 
and calculate the sub-leading corrections to the Wightman function. 



3.2. Higher order corrections to the Wightman function 

By introducing z = —krj and 9(z) = z — (it/A) (1 + 2u), the Bessel and Neumann functions 
are asymptotically expanded around z = 00 (i.e., 770 = —00) in the form 



N v (z) 

Ju+l(z) 



N u+1 {z) 



— < cos 9(z) 

71 Z 



— < sin 9(z) 

71 Z 



— < sin 6(z) 

71 Z 



— sin viz) 



8z 2 
8z 2 



sin 0(z) 



+ cos 



cos 



2 1 



2z 
,2 1 



7T Z 



cos0(z) 



+ cos 9(z) 



^ ^ + cos6»(z) ?- 

2^ z 

x i- 2 -\){- 2 -l) 



2z 
,2 1 



2z 



+ 0(z~ 



+ 0{z 



V 



2 1 



— sin viz) 



2z 



(v 2 -!)(-+ 1) 



22 2 



+ sin 9(z) 



v 



+ 0{z 



(3.28) 
(3.29) 



(3.30) 



(3.31) 
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Then, w has the expansion 
'd-l 



Uq = 



■voY 



'2k 



+ u + iu{r) ) rjo Jy{z) - z J 1+V {z) 



7T 



! + j ¥ - (" 2 - l) (" 2 - l) - i - ( d ~ 2 ) " 2 ^ 2 ) 



2z 



d -2 (v 2 - |) (d - 4) - 2m 2 



(3.32) 



4z 8z 2 

v can be obtained from this expression by replacing 9(z) in (I3.32p with 9(z) — n/2 (see 

(^M-<^MY 

d-l 



V = -(-Tjo) 2 



+ u + iu(rj ) r) N v (z) - z N 1+V (z) 




l *?-(^-\) {v 2 -\)(v 2 -\-{d-2)-2m*Y 
1 2z 8z 2 

d-2 (v 2 -\)(d-A)- 2m 2 



— i- 



Az 8z 2 
With the u and Vq given above, the wave function (13. 23f) can be calculated as 



0(z- 6 )\. (3.33) 



x < e 



-i9(z) 



1 + i ¥-(- 2 -i) (- 2 -l)(- 2 -!-(^-2)) 



2z 



iz' 



d -2 {v 2 - I) (d - 4) - 2m 2 



Hi 2 \-k V ) + 0(z- 3 )\, (3.34) 



from which the Wightman function for each mode is obtained as 



+ [H^i-k r,) Hl 2 \-k rf) + Hl 2 \-k V ) H$-\-k rf)] 



d-2\2 1 



+ H^{-k7 1 )H^\-k7 1 l ) e" 2i 



+ Hl 2 \-kr ] )Hj 2 \-k V ')e 2ie ^ 



.d-2 , (v 2 -l)(d-3)-(^Y 



m 



Az Az 2 
.d-2 , (^ 2 -l)(rf-3)-(^) : 



42 



4z 2 



0(z~ 3 
(3.35) 



Here, the last two terms include a rapidly oscillating factor e ±2ld ^ oc e T2lfcr?0 , so that they 
can be dropped from the expression because they do not contribute to the integral ( I3.26P 
for a sufficiently large \rj \ . 
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We thus find that the Wightman function has the expansion of the form 

oo 

G + (x, 770) = E £ c $ G %>( x > ^ Vo) . (3-36) 

n=0 <r=±l cr'=±l 

where 

G (n) 7T [(-q) (-,/)] ^ 

4 (27r)— (-?7o) n | a; - x'| 2 

x f dkk^~ n Ji^{k\x-x'\)H^\-krf)Hi^~\-kr]'), (3.37) 
io 2 
JO) _i (0) _ (o) _ n (i) _ n 

~~ 1 ' -1,-1 ~~ °1,±1 ~~ U ' C <T,cr' ~~ U ' 

c%=4 2 » 1= (^) 2 , cg = c< 2 !,_ 1 = 0. (3.38) 

Note that the n = 1 terms totally disappear as a consequence of ignoring the highly oscil- 
lating terms. As shown in Appendix [Bj the integral (I3.37P can be performed analytically for 
arbitrary n . In the next subsection, we consider the expansion to the order n — 2, which 
is the first nonvanishing sub-leading contribution. We then only need to take care of terms 
with a' = — a since terms with a' = a do not contribute to the Wightman function at least 
for n < 2 (see (ETSgjl ). 

We close this subsection with a comment that the state |0^ ) must break de Sitter invari- 
ance since it introduces an extra scale rjo as a cutoff. This can be seen from the fact that the 
n th -order corrections to the Wightman function, G^l^x, x') , are not de Sitter invariant for 

n 

n > . In fact, the factor [(— rj) (— ?/)] 2 / (— r] ) n in (IB.7j) breaks de Sitter invariance although 
u a - a is invariant. 



3.3. Calculation of T 

According to the expansion (13. 36ft of the Wightman function, the derivative of the response 
function, (I3.6p . has the expansion 



Wr,r i; r ) =5>- n ^°) ^ n \AE,r,n) 

n=0 

'd- 2\2 



t%(AE, r, n) + e" 2 (— ) (^) Yl HUM, r, n) 



cr=±l 



C(e- 4(r - To) ), (3.39) 
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where r = — log(— ?7o) and 

cr=±l 

d S e- iASs G^ (T (x(r + S ),x(r);r /0 ) 

(r-ri) 
(r-ri) 



(3.40) 



3 n (r-r ) 



+ e" 



dse- l *»'G%LMT),x(r ~ *);»*>) ■ (3.41) 



The integration can be analytically performed as is shown in Appendix [Cj and we obtain 



' r a, — (T 



1(7 



2 3+n 7r ^ r (d f l) 



d-3-n -2 v-n-iAE rf _ 3 _ 



2 ' 



-u-n-iAE 



icrG 44 | e 



0, v, 

d-3-n d-3-n , ^-f-n-iAE d-4-n _, 



2 ' 



0, ZA 



-4— n 



2 

d-1 



-iAB 



_ 3 _ n ^-f-n+iAE 



-3— n 



2 ' 



0, v, 



2 ' 2 

ti-v-n+iAE 



d-3-n d-3-n , _ Y-v-ri+iAE rf-4-n , 



2 ' 



0, u, 



d—4—n 



; -(r- Tl )(^-,-n-iAE) ^ £,2,2 | i w e -2(i— n) 



2 ' 2 

^-i/-n+iA£ 
2 

d-3-n i^—v-n-\AE d-3-n 
2 ' 2 ' 2 

-K-n-iAB 



+ 1/ 



-3— n ci— 3— n 



2 ' 



0,1/,-^. 2 

d _ A _ n 



0, I/, 



d— 4— n 



/A 



-K-n-iAB 



-3— n 



-i/-n+iAE d-S-n 



e -(r-n)(^i-i/-n+iA£) ^ _^2,2 | g0 - i7r g-2(T-n) 



2 ' 



0, !/, 



— i/-n+iAE 



d-3-n ri -3-n , -^—v-n+\AE d-4-n , 



2 ' 



0, iA 



d— 4— n 



2 ' 2 

%l-i/-ji+iA.E 



(3.42) 



where G™^ is Meijer's G-function (see Appendix |A|). As is also shown in Appendix [Cj this 
has another expression in terms of the generalized hypergeometric function: 

fl%(AE,T,T i; 770) 



e <jK(AE+if) y 



^=i+v-n+iAE\ „ /^i-i/-n+iAS\ p/^+iz-iAE 



-v-iAE 



2 3 + n 7r ^r(^_l) 
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23+« 7 r i ¥r(^) 



d—l — n 



d-l-n ^+u-n-iAE d _ t _ 



-am ~ ; - -amv / 2 ■ u-x-n , 

e 2 e ^/2' 2 '2^^0-10 



Sill 7T V) 



d—l — n 



3^2 I d±3 + „_ n _ iAj E ; e ° 

1 + V, — 5 



d-l-n ^-v-n-\AE d _i. 



_ — (T 17T s „(T17rl/ / — ^ — y 

e 2 e / 2 ' 2 '2 ^ 

Sin(— 7TZ/J \ \ — V 2 



,<r iO 



I 2 e pf 2 ' 2 '2 CTi0 

sm(7r i/) V 1 + v 2 +» 



d—l — n 



d-l-n ^r-—v—iAE d-l-n 



a — a i7r — = — „cr i7r f / — ^ — y 

e ae t? i 2 1 2 '2 ^ 

Sin(— \ 1 — 1/ -2 



a- iO 



■ d—l — n / d — 1 

-<T 17T 



2 

d-1 



d-l-n ^+u-n-iAE d-l-n 



emu p -{^+u-n~iAE){r-T 1 ) _ / SziZ", 2 , — , + V 

1 + 1/, - 



e 2 e e ^ ^ ^ p 2 ' 2 ' 2 ' • P -2(r-n) 

+ 3^2 | . m+u-n-iAE ' e 







sin(7r v) 


e -ai7r 


d-1- 
2 


± e *i-K V e -(4fL-u-n-iAE) (r-n) 






sin(— 7T i/) 


e oi. d 


2 


e <T ijri/ e -(^+^-n+iA_E) (r-n) 






sm(7r V) 


■ d 


2 


e -cr ivri/ e -(fci-i^-n+iAB) (r-n) 



d-l-n V--^-n-iA£ d-l-n _ 

, pi 2 ' 2 '2 47 - 2 (r-n) 

i _ ^ M^ ra -iAE ) e 

d-l-n ^+v-n+iAE d -l-n , 
I 2 ' 2 '2 -2(r-n) 

+ — * 3^2 I <*±3 +I/ _ n+iA i5 i e 

l + i/, -a ^ 

d-l-n ^-v-n+iAE d_l_ n _^ 



2 ' 2 '2 - . -2(r-n) 

H W \ 3-^2 J±3_ v _ n+i AB i e 

Sin(— 7TZ/) \ 1 — V, 2 " n+im 

(3.43) 

Here, 3F2 is defined from the generalized hypergeometric function 3 F 2 (see Appendix \K§ as 
9 /a 1; a 2 , a 3 _ x = I>i) r(a 2 ) r(a 3 ) /ai, a 2 , a s . \ . . 

3 2 v fti.fta 'V- r(&or(& 2 ) 3 2 \ b u b 2 >*)■ l ^ 44J 



For the case of n = (then we only need to take a = —1), we obtain 



+(0) 
•'-l.l 



e -7T AE y I ^I±^±1^E\ y( ^- u+iAE ^ T,/ £ ^ i +^iAi?^ t,/^- u-iAE" 



. . d±3 +i ,-iA 

1 + z/, -^— 2 — 

/d-1 i^—y-iAE d _i 



, : p 2 > 2 > 2 ■ _ _-2(t-ti) 

8 7T 2 1 (^ j sm(7r v) \ 1 + u, - 



i7 r,g-(^-,-iA E)( r-n) / 2=± -2— ^ , ^ - y 

^ : — d-i _/j ^ — : 3-F2 1 d+i- u -iAE ; e 



8tt— r(^fi) sin(-7rv) \ 1-1/, 



-1 -^^-(^^(r-n) ^ Z£±^ ; d^l +I/ 

PI ■ — 2(r-ri) 

3^2 I M+,+iAE ' e 



8tt— r(^=i) sin(vrz/) \ 1 + z/, 



2 
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I 2 e e 2 1__ 9 2 ' 2 '2 " _2(r-ri) 

H 5=i 7T-T ! ! 3^2 d±3_ y+ i A i? ) e 



8 7r^r(^i) sin(-Tr^) \ i-^: 



2 



(3.45) 



If we here set d — 4 and z/ = 1/2 in (I3.45p . the result obtained in section IV-D of [H] is 
correctly reproduced. Moreover, in the limit r — T\ — > oo, ( I3.45P gives the known result 
which is applicable to arbitrary dimensionality [T4"t [T5~|: 

e -rfABp/'3^+^£\ p/ ^- y +UAg \ T / *^+v-itAE \ r / ^-u-i£AE \ 

f c %ae) = ^ — - — LA — 2 — — LA — i — LA — i — L 

87r^r(^i) 

(3.46) 



where we have restored the curvature radius I temporarily. Note that this function satisfies 

e- 2 * eAE . (3.47) 



P«(AE) __ 27v£AE 



F°i(-AE) 

For later convenience, we here give the large \AE\ limit of the function 

A(AE) = F cq {AE) + F cq {-AE) . (3.48) 
This can be calculated by using the asymptotic form of the Gamma function (8.328-1 of 



lim \F(x+ iy)\ ~ V2^e~^ lvl \y\ x s (x, j/GE), (3.49) 

|y|-Kx> 



and we obtain 

id— 3 



I A .El 

A(A£) ~ - ' J n/ . n ( |A£| -+ oo) . (3.50) 



4. Thermalization of a two-level detector in de Sitter 
space 

4.1. Two-level detector 

We first summarize the assumptions we have set in deriving equations fl3.4p -( l3T7|) : 
1. The scalar field 0(7/ , x) is in the ground state at an early time rjo . 
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2. The detector starts interacting with the field at a later time rji > rjo), when the 
total density matrix takes the form 

p tot (n) = p(n)® K)(0 VO \. (4.1) 

3. jo(r') in the integral of (I2.36P can be replaced by p(r) due to the singular behavior of 
the Wightman function around r' = r . 

In this section, we consider the case where the detector is a two-level system with en- 
ergy eigenvalues E\ and E 2 with AE = E2 — E\ > . We make the following additional 
assumptions: 

4. The energy difference is much larger than the natural energy scale of de Sitter space, 

ae > r 1 = 1 . 

5. The initial distribution of the detector, p{ri), is averaged over the initial proper time 
Ti for the duration At > 1/ AE in order to describe distributions at different energy 
levels to a good accuracy. As a consequence, the off-diagonal elements of piji) can be 
effectively set to zero because they are oscillatory during the time. 

6. The operator /1 is off-diagonal with respect to the basis \m) (m = 1,2), 

A*=«mHn»=(° ^ 



7. The detector is an ideal detector, in the sense that it gets adjusted to its environ- 
ment very quickly. We see below that this condition is expressed as the inequality 
\ 2 \^ 12 \ 2 {AE) d ~H^> 1. 



4.2. Relaxation to the Gibbs distribution 

Under assumption |6] for a two- level system, one can easily show that ( I2.37P decomposes into 
the diagonal and off-diagonal parts. The off-diagonal part is given by the linear differential 
equations 

iAE - |// 12 |' 4 J*(0; r, n; rj ) lAi^i^ T i Til ^0) 

/ill j"(0; r, n; 770) - iA-E - \fi 12 \ 2 £(0; r, n; 770) 

MM,) (43) 
\P2i(r)J 
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which are integrated to 



^ T »)=Texp(/V.M(r',, 1 )) f^' 

yp2l{T) V ri 7 \P2l(Tl). 



(4.4) 



Since we can effectively set pui^i) = P2i( r i) = due to assumption [51 we can also set 
Pi2(t) = P2i(t) = (assuming that Ai(r',Ti) changes slowly when taking the average over 

Tl). 

As for the diagonal part, using (I3.42p together with assumptions |4] and |5l one can easily 
show that those terms in JF which explicitly depend on t\ disappear after taking an average 
over T\ , and T comes to take the form 

pri+Ar 

Ai 



1 r +A7 " dr , -p^, T) r ,. = p^ AE) + e -2(r-r ) ^(2) (AE) + Q^- 

' J Tl 



4(r— r )> 



(4.5) 



where ^^(AE) is given by (13.461) and 
r(2) ,.„, /^-2x2 1 



J 7 ^ j (A_E) 



4 y 2 3+«7r^r(^ i ) 



X 



-3— n 



-u-n-iAE 



-3—n 



iGgle- 



-iG^le— 



2 ' 



0, !/, 

rf— 3— n 
2 ' 



2 ' 2 1 

<*zi_ l/ _ n _iA.E 
2 

^-u-n-iAE d _ 3 _ n 



-v-n-iAE 



+ (AE —AE) 



0, 1/, 

d-3-n _d-3-n i ^-y-n-iAE d-4-i 



2 ' 2 
0, !/, - 



-4-n , .. V-^-n-iAg 



d— 3— ri d— 3— n 



/A 



2 ' 2 
0, V, 



2 ' 2 

^-u-n-\AE 
+ V 

1 3 O 



(4.6) 



We thus find that the master equation for the diagonal elements takes the form 



Pn[T) 
,P2 2 (t) 



-w + (t;t ) w_(r;r ) \ f p n (r) 
w+(r;r ) -w_(r;r )/ \P22(t) 



(4.7) 



where 



w±(t;t ) = A 2 |/i 12 | 2 [P\±AE) + e- 2 ^ T - To) ^ 2 \±AE)] 
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From (I3.47P we see that, as r becomes large, the density distribution p mm (r) approaches the 
Gibbs distribution at temperature l/2n£ (see ( I2.43P ). 

-2n£E m 

P2 m = -^- (Z = e-^+e-^). (4.9) 
In order to investigate how the detector relaxes to this equilibrium, we expand p mm (r) 

as 

p mm (r) = p^ m + Ap mm (r) (4.10) 

and keep small quantities to the first order, assuming Ap mm /p^ m <C 1 and e~ 2<r ~ ro ) <C 1 . 
Then (14 .7p becomes 

^Ap n (T) = A 2 |/i 12 | 2 [- e ~ 2 ^o) j&)( AE ) p% - ^(AE) A Pll (r) 
dr 



+ e -^ T - To) P 2 \-AE) p c £ + P^(-AE) Ap 22 (r)] . (4.11) 



Noting that P^ 2 \AE) is an even function of AE (see (14.61) ) and that Apu = — Ap 22 , we 
can rewrite (14. lip to the form 

-^Apn(r) = -A 2 |pi 2 | 2 \e~ 2(T ~ To) P 2 \AE) tanh(7rA£) 
dr 

+ (F*(AE) + j^(-AE)) Apu(r) ] . (4.12) 
The general solution to this equation is given by 
Ap n (r) = e- A2 l^l 2A ( AE )(^)Ap 11 (r 2 ) 

tanh(7rA£)( e - 2 ^ - e'* l/^^K^Mto-^)), 



A 2 |p 12 | 2 A(A£) -2 

(4.13) 

where A(AE) = J^^AE) + J rcq (— AE) (see (13.48}) ) and r 2 is an arbitrary proper time after 
T\ (r 2 > ri). We see that Ap(r) certainly approaches zero as r increases. 

Moreover, since A(A£) ~ | AE\ d ~ 3 as A_E £ > 1 (see (I53D|I ) and A 2 |pi 2 | 2 (A£) d " 3 £ > 1 
as stated in assumption 0, the damping terms proportional to e - ^'^ 12 ' A ( AE ) ( T_T 2) rapidly 
disappear from the expression, leaving the term proportional to e~ 2 ' T ~ T2 ' ) , 

Apn(r) ~ - tanh(vrME) e -^)/' , (4.14) 

where we have restored the curvature radius Note that the coefficient ^^(AE) / 'A(AE) 
does not depend on details of the detector (such as Api 2 ). Since only those rapidly disap- 
pearing terms depend on such details, one may say that the detector under consideration is 



22 



an ideal detector, in the sense that it quickly looses its own nonequilibrium properties and 
gets adjusted to its environment. Then, the remaining damping term in (j4.14p is interpreted 
as representing the relaxation process of the media surrounding the detector to the equilib- 
rium at temperature l/2ir£ with the universal relaxation time £/2. Note that the reason 
why the relaxation time is given by £/2 (not by £) is the vanishing of the n = 1 terms in 
( I3.36p . This means that the relaxation time is still given by £/2 for any ideal detector which 
is not necessarily a two-level system. 

5. Conclusion and discussions 

In this paper, we have considered an Unruh-DeWitt detector staying in the Poincare patch 
of de Sitter space. The main difference of our setup from that given in the literature is 
in that the scalar field (weakly interacting with the detector) is not in the Bunch-Davies 
vacuum (nor in the a- vacuum), but in the instantaneous ground state \0 Vo ) at finite time 
rjo . In order to deal with such circumstances, we first derived the master equation which 
describes a finite time evolution of the density matrix of an Unruh-DeWitt detector in 
arbitrary geometry. We then applied the framework to de Sitter space and showed that 
there always exists a relaxation mode with the universal relaxation time £/2 , irrespective of 
details of the detector. 

We here should stress again that we are not considering just the thermalization process 
of a detector dipped in a thermal bath. In fact, the state \0 Vo ) is different from the Bunch- 
Davies vacuum, and thus, the detector initially should behave as if it is in a medium which 
is not in thermal equilibrium. As time goes on, the detector comes to behave as if it is in a 
thermal bath, since the difference of the state |0 W ) from the Bunch-Davies vacuum becomes 
irrelevant at later times. 

In reality, there can be many relaxation processes for such a detector, some of which are 
simply the processes where the detector becomes adjusted to its environment. However, these 
processes usually depend on details of the detector, and thus can be neglected by assuming 
that it is an ideal detector which quickly responds to changes in its environment. The terms 
proportional to e~ x2 ' Ml2 A ( AE "> i T - T ^) in (14.131) actually represent such processes. On the other 
hand, the damping terms proportional to e~ 2( - T ~ T2 ^ e in (I4.14p necessarily survive even for an 
ideal detector, and thus describe the nonequilibrium dynamics of the surrounding medium 
relaxing to the thermal equilibrium corresponding to the Bunch-Davies vacuum. We thus 
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can conclude that the relaxation time £/2 is a quantity representing the nonequilibrium 
properties intrinsic to de Sitter space. 

In the previous studies, thermal properties of de Sitter space have been investigated 
only in the context of equilibrium thermodynamics, and its nonequilibrium properties have 
not been studied well. The relaxation time calculated in this paper may serve as the first 
example related to the nonequilibrium dynamics intrinsic to de Sitter space. 
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A. Meijer's G-function and generalized hypergeometric 
function 

Meijer's G-function G™ q n (0 < m < q, < n < p) is defined by (see 9.302 of [IB] for details 
on the choice of the contour C) 

a t ,...,a p ^ _ f ds UT=i T ( b J- s ) ITU r (! ~ «i + s ) ^ ^ 



P ^ h : ...,b q J J c 27ri Y[] =m+1 r(i - b, + s) uU+i r (% - s ) 

This function is invariant under arbitrary permutation of a set {ai, . . . , a n }, {a n+ i, . . . , a p }, 
{£>i,..., b m }, or {b m+ i, . . . , b q }. As a special case of the G-function, the generalized hyper- 
geometric function is defined by (9.34-8 of [i~6] ) 



Q!,..., Op \ _ Ui=l T ( b i) n \ 



p g[ bu-.-A )~nun*i) p ' q+1 



1 — Ci\ , . . . , 1 — dp 

0,1 1-6, 

nLir(6i) f ds rHiE =1 r(o i + fl ) 



nr =1 r(a 4 ) j c 2ni uum+s 



-z) s . (A.2) 
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For convenience, we define the following functions using the generalized hypergeometric 
functions: 



Oi, . . . , a % 



, . . . , Op ^ 

, • • • , bg ' 



(A.3) 



If no two bi (1 < i < m) differ by an integer, Meijer's G-function G^Y 1 with p < q, or 
p = q and m + n > p, or p = q and m + n = p and |z| < 1, can be expanded by using the 
generalized hypergeometric functions [17], 



p,q 



, . . . , dp 
,...,b q 



m n m 



E 



l-ai + b k , 



1 — a„ + bi. 



bk, ...,*,.••, 1 — b n 



(-iy- m - n z) , (A.4) 



where * means that the k th term has been omitted. In particular, G 3 ' 3 is expanded as 



ai,a 2 ,a 3 s 

bl T(b 2 - 6Q r(l - 6 2 + 6i) p /1-ax + Oi, 1 - a 2 + 6 X , 1 - a 3 + &i . _ z 
r(a 3 - 6i) T(l - a 3 + 6i) 3 2 V I-62 + 61, I-63 + 61 



fe2 r(6 1 -6 2 )r(i-6 1 + 6 2 ) ~ 

T(a 3 - 6 2 ) T(l - a 3 + 6 2 ) 3 2 
Furthermore, we can show 



1 - ai + 6 2 , 1 - a 2 + 6 2 , 1 - a 3 + 6 2 
1 - &! + 6 2 , 1 - 6 3 + 6 2 



(A.5) 



G 4 ' 4 ( z 



a li a 2, a 3i a 4' 

6i,6 2 ,6 3 ,a 4 > 



6l r(6 2 -6 1 )r(l-6 2 + 6 1 ) £ 



+ Z 



T(a 4 — 61) T(l — a 4 + bi) 

r(6 1 -6 2 )r(i-6 1 + 6 2 ) 



62 



T(a 4 - 6 2 ) T(l - a 4 + 6 2 



4-^3 



bi r(6 2 -6 1 )r(i-6 2 + 6 1 ) 

T(a 4 — 61) r(i — a 4 + 61) 

b2 r(6 1 -6 2 )r(i-6 1 + 6 2 ) 



+ z ( 



1 - ai + 61, 1 - a 2 + 61, 1 - a 3 + 61, 1 - a 4 + 61 _ 
1 - 6 2 + 61, 1 - 63 + 61, 1 - a 4 + 61 
1 - ai + 6 2 , 1 - a 2 + 6 2 , 1 — a 3 + 6 2 , 1 - a 4 + 6 2 
1 - 61 + 6 2 , 1 - 63 + 6 2 , 1 - a 4 + 6 2 
1 - ai + 61, 1 — a 2 + 61, 1 - a 3 + 61 
I-62 + 61, l-6 3 + 6! 
1 - ai + 6 2 , 1 - a 2 + 6 2 , 1 - a 3 + 6 2 



T(a 4 - 6 2 ) r(i - a 4 + 6 2 ) 3 2 V 1 - bi + 6 2 , 1-63 + 62 
where the following identity (which can be shown from ( 1A.2I) ) has been used: 

'ai, a 2 , 03 



.ax, a 2 , a 3 , a 4 

61,62,04 



3^2 



61,62 



(A.6) 



(A.7) 
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From the above identities, we can show the following relations for a = ±1 : 



_ n _ 3 n- i^+u+iAE _ d _ n _ - 



2 ' 



o, 



e- ,if Y sin( 7 r^ Ii ) 



d-l-n d-l-n , ,, VH-^-n-iAg 



Sin 7T 1/ 



3-^2 



2 ' 2 
1 + V, 



d+S 



sin 



/ d-l-n-2iA 
^ 2 ) 



d-l 

d—l—n d—l—n , 2 



+u-n-iAE 



u-n-iAE 



e ai0 x 



sin — 7T 1/ 



3-^2 



±tl- u - n -iAE 
1 — I/, » 



n cr iO 



4 | 6 X 



d _ n _ 3 n- Y+^iAE _ d-n-3 

o 1 Is rt 1 Is 



-71 — 4 ' 



2 ' 



u, 0, 

-I-n\ 



2 ' 

n-^+u+iAE _ d _ n _ 4 

O 1 ^ O 



Sin 7T ^ 



- — 3F2 



d -l-n d-l-n 1 — +^n-iAE 

2 ' 2 ,. ' 2 .^io T 

1 + V , — n 



COs(tT <*-l-n-2v \ ^ 
+ \ ( \ 1 3F2 

sm(— 7T v) 
from which we obtain the formula 



d—l—n d—1- 



2 ' 2 



-v-n-iAE 



1-U, 



d+3 



v-n-iAE 



e CTi0 a; 



d— n— 3 
2 ' 



-n-3 «-¥+iAE+i'' 



i/, 0, 



2 ' 



2 
1 



_ ■ _ d — n 



d-n-3 _ d-n-3 n-^+u+iAE _ d-n- 4 ' 



2 ' 
n-t^+HiAE 

1/, 0, 2 — 9 



hu- 



-n— 4 



^i+y-n-iA_B 



sm 7T 1/ 



d —l—n d—l—n 1 2 

2 > 2 ' 2 

1 + ^ — 5 



e ai0 x 



d—l—n d—l—n 



- V, 



-v-n-\AE 



+ 



sin — Ti v) 



3 r 2 



2 ' 2 > 2 

^-v-n-iAE 

y — v, — „ 



; e aio x 



Another useful formula can be obtained by using the formulas [T8 

3-^2 ( , , 5 2 

V 61, 62 

F(ai) T(a 2 - ai) T(a 3 - a x ) /ai, ai - 61 + 1, ai - b 2 + 1 

,—2) zrrz ^ — ; 3^2 I - . \ Z 



r(6i -ai)r(6 2 - ai) 



a\ — a 2 + 1, Q-i — a 3 + 1 



_l_ (-^)" a2 r ( Q 2) r ( Q i ~ Q 2) r ( Q 3 ~ gg) /a 2 , a 2 - h + 1, a 2 - b 2 + 1 ^ 
T{bi- a 2 )T(b 2 - a 2 ) V a 2 — ai + 1, a 2 — a 3 + 1 



r(a 3 ) r(eti - a 3 ) T(a 2 - a 3 ) /a 3 , a 3 - 61 + 1, a 3 - 6 2 + 1 



r(6i - a 3 ) T(6 2 - a 3 ) 



3^2 



03 — ai + 1, a 3 — a 2 + 1 



; z 



a\ — a 2 Z, ai — a 3 ^ Z, a 2 — a 3 ^ Z, and 2 ^ (0, 1) 
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from which we can show 



2 ' 2 '2 

1 ±u, -a ^ 



_ 2j 2V 2 ^ j sinl7r— 



V 2 

X 



7T 

d-l-n d-l-n-,-,, ^T^-iAS 



■ /■_ d-l-n T 2i^ / d-l-n d-l-n -■- 

+ 2 3T7 \ 3^2 | d+3 Tl/ - iA E i Z 



sin(^7rz/) I i T ^ 



-l-n d-l-n i V^-iAS 



2 



„• ^ d-l-n A / a-i-n a-i-n i 

d-l-Ti±ai/ Sin(7T j J 9 [ 2 ' 2 ' 2 

sin(±7ri/) \ i ± „ _a_ 



Using this equality, the following formula can be shown for x > 0: 



3^2 d±3 +I/ _ n+iA £; ' e X 



d-l 
2 



su^tt i/; V 1 + v. 



2 



' X 2 V 2 ' ~ / 2 ' 2 ^' 2 _ CT i0 



3^2 «- v - n +iAE 5 e X 

riwt) r( V^ 1 " ) r( ^-T +IAE ) r(g±f^) r(fe_^) 

7T 

1 / d-l i ,, : A rp\ _ ; _ c d-l-n , , \ /d—l—n d—l—n i ^^-+v— iAE 

H r-7 s 3-^2 d+3 +iy _iAE ' e X 

sm(7rz/) y i + „ 9 2 + ^ 

.If <*=!_„_< A Jtf -„ W*- 1 -" -,/) /d-l-n d-l-n ifi-u-iAE 



i / w— i ,, ; a p\ j /• a—i—n \ / a—i—n a—i—n 

( — -u-iAE) e -a^{— „) / __ _ ^ _ i 

H r-T s 3-^2 d+l-^-iAE ? e x 

Slll( — 7TU) \ i _ y 



(A.14) 



B. Wight man function 



In this appendix, we show an analytic expression for the integral given in f !3.37p . By using 

. n ( 3+tT ) 
the identity K v (e lcr 2 £) = — (icr7r/2) H y 2 (z) with a = ±1 , we first obtain the relation 



dkk^ +1 ~ n J^s(k \x - x'\)H^\-kri)HF* \-kri') 



'o 
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A I POO 

e (-+-') ^ / dk k ^+i-n \ x - x'\) K u (-e [ ° i fc 77) K,(-e iCT ' ? fc 77') 

^ JO 2 



(B.l) 



Then, by using the formula 

poo 

/ dk k" +l - n K u {a k) K u (b k) J At (c k) 



2^— 1— n c /i 



r(/i + 1) 6 2 « 

x (a/6)" r(0 r(c + 1/) r(-v) f 4 (£, £ + 1 + 1 + ^; -c 2 /6 2 ; «7& 2 ) 
+ (a/6)-" r(0 r(e - 1/) i» f 4 (£, £ - ^; 1 + n, 1 - f; -c 2 /6 2 ; « 2 /^ 2 ) 

[e = /i+l-(n/2), i/^Z]. 
we obtain the following expression for G ™l,(x, x'; 770) : 



(B.2) 



G^],(x,x';r/o) 



2 2+n 7T a e 17Ta — 



r(^) v-77 



x F/ 



77'/ 2 

d — 1 — n d — 1 — n d — 1 



e iCT '° I a; - sc'|" / e ian (-rj)' 



' 2 ' 



1 + 1/: 



-77') 2 ' e if7 '^(-r/) 2 



+ (!/->■ -i/J 



(B.3) 



Here, F 4 (a, 6; c, d; x, y) is Appell's hypergeometric function given by 



00 00 



F 4 (a,6; c,d;x,y) = 



m=0 n=0 

(x) n = T(x + n)/T(x) 



(^Om+n (6)m+n 

m\n\ (c) m (d), 



x' 



y n [N 1/2 + M 1/2 <i], (B.4) 



(B.5) 



For Z(x, x') = [r] 2 + r]' 2 — (x — a?') 2 ] / (27777') < 1, the series in (1B.4[) does not converge and the 
analytic continuation should be performed (see 9.185 of [TB] for the integral representation 
of the Mellin-Barnes type). 

In the limit \x — x'\ — > 0, one can show, by using the identity 



F 4 (a,6; c,d;0,y) = ^ 



n=0 



i a )n (6)n n „ fa, b 



(B.6) 
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that the Wightman function takes the form 



lim G„l,(x,x r ) 

\x-x'+0 ' 



-7i y e^'l^r ^f^^ + u,^^ e ia7r (-rj) 2 

2 F\ 



rj'J sin(7rz/) \ 1 + v ' e la ' n (— rj') 2 
-riy e-^'y^r - - ^, . e i<J7r (-r/) 2 



+ — 7 — -7 r 2F1 



-77V sin(— ix v) \ 1 — v ' e lcr ' w (—r]') 2 

(47r)fr(^i) (-»»)» 2 \ ; 2 

(B.7) 

with 

U ™' ~ 2e i (-+-')(— o) (_^/) • ^-») 
Here, in the last equality, we have used the following formulas (see 9.132-2 and 9.134-3 of 

HE): 

2Fl \ 7 > Z ) ~ TWTb - *)^ Z) " Fl \ a + l-P '* 



+ 



r( 7 )r(q-/3) p ,0,13 + 1-7 -i\ 

I»r( 7 -/3) 2Fl l /3 + 1-a '* J 

[|arg2:|<7r, a - /3 7^ ±m , m = 0, 1, 2, • • • ] , (B.9) 

In order to understand the behavior of the Wightman functional 
G + (x(t),x(t - s); ?7o) 



10 The Wightman function has the following asymptotic form for large s : 



G + (x(t),x{t - s); 7? ) r. 

4 7T 2 



d + l 

Ait— 



which becomes oscillatory for the heavy mass case m > (d— l)/2 . The behavior for small s is given by 



d-2 

e 



G+(*(t),x(t - s); Vo ) ~ / 2 > s~^) 

Air 2 
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(7=±1 



C^ (T (x(r),x(r - s); r/ ) 



(B.ll) 



we give plots of the Wightman function for a few typical cases of d = 4 with z/ = 20 i 
(heavy mass case), 1/ = 0, and v — 1.4 (light mass case). It is clear from these examples 
that, for the cases other than the light mass case (i.e., for the cases described in figures [3] 
andH|), the Wightman function G + (x(r),x(r — s); Tjo) takes significant values only around 
the coincident point s = 0. On the other hand, as can be seen in figure El the Wightman 
function has a correlation with a longer range as the mass m decreases. Thus, in the light 
mass case, we may not be able to neglect memory effects and the Markovian approximation 
used in (I2.37P may not be valid. 




0.5 1.0 1.5 2.0 2.5 3.0 




Figure 3: The real part (left) and the imaginary part (right) of the Wightman 
function G + (x(r),x(r — s); 770) for d — 4, v — 20i (heavy mass case), and 
r - To = 50. 



0.2 0.4 0.6 0.H 1.0 



0.2 0.4 0.6 0.8 1.0 



Figure 4: The real part (left) and the imaginary part (right) of the Wightman 
function G + {x{t) ) x{t — s); t] ) for d = 4, v = 0, and r — r = 50. 
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Figure 5: The real part (left) and the imaginary part (right) of the Wightman 
function G + (x(r), x(t — s); 770) f° r d — 4, v = 1.4 (light mass case), and r — r = 
50. 



C. Derivation of $2321) and ( ESp 



In order to calculate P^L a for the trajectory (I3.2p . we should consider the limit |a; — x'\ — > 
in the Wightman function G + (x, a/). In this limit, the n th -order Wightman function ( 13 . 3 7|) 
takes the form 



< 3-<r 1 



2 (47r) ^r(^i) Jo 

(C.l) 

Then, the first term in (13.41)) is rewritten to 

= n(T " T(,) / dse- iAEa G^L a (x{T + 8),x(r)) 

J-It-tl) 



e 



2 M ¥r(ti) 7_ (T _ n) 



/ ds f^(x/ y )^x^y¥-»/ff^x)/ff?>(y) 

J-(t-ti) JO k 



/»oo 2_ /*e( T r i) y 



2 (4vT)^r(^i) Jo ^Jy 

| - ' / ix&«» H?*Hx) / dyy^-»— H^\y). (C.2) 

2 (47r)— r(^-) Jo Je-tr-' i)x 

Here, in the first equality, we have defined new variables x = k e~ ( - T+s ' ) and y = k e~ T , and in 
the third equality, the order of integration has been changed. In a similar way, the second 
term in (13.411) is rewritten to 

,n(r-T ) f ds e" 1 AS * G^L (^(t) , 2 (t - S ) ) 



7t i r~ Ti . r°° dk 



2 (47r) ^r(^i) Jo 



/ T " T1 d. r ^ (x/y)— y--x- irf^(y) irf^(x) 
Jo Jo fc 
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7T 



2 (4vr)^r(^i) 
1 



7T 



2 (4vr)Vr(^i) 



I q_|_ 



I o 



.0- a) 



-(t- n ) x 



dyy^+^tff^y), (C.3) 



where we have defined x = ke ( T ^ and y = ke r , and again the order of integration has 
been changed in the third equality. 

The y integration in (1C.2P and (IC.3P can be performed as follows: 



dyy^-^^tff^y) 



2i<J _ iTTK 

7T 



=F — e + 2 



fox 



dyy^-"T iA ^(e^y) 



7T 



1(7 

=F — e ~ 

7T 



(ax) 2 /-(bx) 



dF 



dFF" 



ax 2 



y — 2 — ^ e^oxvy - (a -►&) 



2"" 1 icx 



=F- 



7T 



/ pTCTlTT / ^2 



*A 0, 



(a -»■ 6) 



In the last equality, we have used the following formula (see 6.592-2 of [16]): 



(C.4) 



2 ~ A 

u, 0, I - A - 1 



2^ 1 a~ 1/ i»r(l-z/) 1 F 2 



A + 1 - f 
l-i/,A + 2- |'T 



-2-"- i ^r(-z/)r(i + z/) 1 F 2 



A + 1 + ^ a s 



1 + v, A + 2 + | ' 4 



Re A > -1 + 



\Reu\ 



(C.5) 



If 6 0, the condition Re A > — 1 + (|Re^|/2) [<^> (d — l)/2 > |Rez/| + n] is not necessary. 
In order to perform the x integral in (1C.2j) and ( 1C.3|) . we use formulas 7.821-1 and 7.821-2 

of eel 



2,1/wr 



&xx a ~ l J p (x) Glj( 



dxx a ~ l N u {x) Gl'l( U 



h, b 2 , h 
a i 



2 a ~ 1 G 2 / 3 l u 



a+v—2 a— v— 2 

2 ' !> 2 

a+v— 2 a— ^— 2 



(C.6) 



2 ' 2 ' 1> 2 

a— v— 3 



&1, &2, &3, 



(C.7) 
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from which we obtain the formula 



h, b 2 , h 



Gjs ( w 



g+i/— 2 a— ^—2 

2 ' !' 2 



=F icr G4 4 ( a; 



2 a— 1/—2 a—v—3 

95 95 1 ' 9 



6l , 62, 63, 



a— v— 3 



(C. 



Using this, we obtain the formulas 

7T 1 



2(4vr)^r(^) 
±ia 



dxx i_3 ±iA ^(^) (x) 



(-r-n) x 



dyy 



# 2 ; (y) 



2 3+n 7r ^r(^i) 



X 



d-3-n ^-v-n^\AE d _ 3 _ 



2 ' 



0, ZA 2— 



d—3—n 



2 

-3— n 



+ V 



-u—n^fiAE d—4- 



=F 



icr G 4 4 I 



2 ' 2 
0, 1/ 



+ v 



d—i—n 1 2 

2 "r" ' 2 



1/— n=FiAi? 



-(r-n) J ^2,2 f i ^ e -2(r-n) 



d— 3— n 



-v-n^\AE d-3-n 



2 ' 2 

0, v, 



^--y-n^\AE 



-3— n d— 3—n 



2 ' 



-i/-n=FiAS d-4-n 



0, U, + !/, 



2 

d-l 



(C.9) 



Finally, by adding ( 1C.9I) with different signs, equation (13. 42ft can be derived. 
By using identity (lA.lOj) . equation (13.421) can also be written as follows: 



T. 



(n) 

(7,— (7 



/d-l-n ^+u-n-iAE d-l-n 

~— fT 17T ^ „ — U 17Z U I — - 

e ae s / 2 ' 2 '2 



+ i/ 



sin 7T i/ 



d — 1 — n 
g — (T 17T j g CT 17r V 



3-^2 



1 + V, 



ri+3 



iO 



+v— n— iAE 



d-l-n ^-u-n-iAE 



sin — 7T 1/ 



e <J i7T d | " gCT i7TI/ 



3-^2 



2 ' 2 
1 - I/, - 



— V 



sm 7T 1/ 



3-^2 



d-l-n V-+^-n.+ iAE 

' 2 . „-ffi0 



2 ' 



1 + v, 



it3+v-ri+iAB 
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3< r iv i=$=» ( d-l-n ^-v-n+lAE d-l-n _ 



3^2 I *fc3_ v _ n + i A£ 5 e " 



2 



ff i ffVp -(JS f i+ l /-„-iAB)(r-7i) ^ Vz^! — + V 



d-l-n ■ ,, rd-l,„ „ ;ae<W~ /d-l-n -5 — h^-n-iAS 



I C Z C C ' ^ ' ' Z ?| 2 ' 2 ' 2 ' " . „eriO .-2(t-tO 

+ — n 3^2 I d+I +1 ,_ n _ iAB ' e e 



+ — — 3^2 I d+3_„_ n _ iA £ > e e 







sin(7r z/) 






d-1- 
2 


V'"e-(T- 


v-n-iAE) (r-n) 






sin(— 7T v) 


3 




-1 — TJ 

2 




-n+iAE) (r-n) 
3^ 






sin(7r 1/) 




■ d 


-1-71 

2 


/ d— 1 

g-<T 17TI/ g-(-2 


I/— re+iA-E) (r-ri) 






sin(— 7T v) 


3 



1 + u. 



2 



-1-n ^-f-n-iAE d_l_ n _ 



2 



1-n V+^-n+iAg rf_i_ n , 

2 ' 2 '2 CTi0 _2(t-ti) 

+ — ^ 3^2 1 4±i +u _ n+iAE , e e 

1 + -a 2 

-1-n ^-^-n+iAE d-l-n _ 
2 ' 2 '2 _ CTi0 2(r-n) 

_ M- y - n+ i AE , e e 

(CIO) 

Furthermore, this can be rewritten to the following form by using flAlij) with x = l: 



e «(AE+f) p/^j^-n+iA£\ p / i|^-n+iAg \ p/ ^+pAg N p / ^-piAg 

= 23+«7T^r(^i) 
1 

+ 2 3+n 7r Vr(^ 1 ) 

_ CTl7r ^ i7r;/ / d-l-n ^+*-n-iAB d-l-n , 



3-^2 1 *±3 +zy _„_iAB 

sm^7rz/j \ \ \ v 2 

v^i^ll inu /d-l-n ^-u-n-iAE d-l-n. _ 

e 2 e p 2 ' 2 ' 2 CTi0 

~t \ 3-^2 4+3 _ iAE , e 

sm(-7rz/) y x-v, 2 2 

I e 2 e P 2 ' 2 '2 ^ o-iO 

H w \ 3-T2 d+3, iAE , e 

sm(7r i/) ^ x + ^ 2 +u 2 

CTi ^^i^ . /d-l-n ^-u-iAE d-l-n _ 

, e 2 e 9 2 ' 2 '2 ".-id 

H r-7 s 3-^2 d+3 iAE , e 

sm(-7rz/) y 1 - z/, 2 2 

d-l-n ,d-l , .. „ !A e.w_ _^ /d-l-n ^+u-n-iAE 



I ^ C C ^ P 2 2 ' 2 . _«riO.-2(T-TL 

H ~< \ 3-^2 iAE ,e e 

sm(7r v) \ 1 + 1/, 2 + 2 

fflff ii=J=2 (Srl-^n-iASJCT-n) / d ~ X - n fci^Z-^g d-l-n - v 

"1 w \ 3-^2 d+3_ iAE , e e 

sm(-7ri/) V 1 _ y 2 ^ " ^ 
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I C C P 2 2 2 . -cr 10 -2(t— 71 

+ sin^ ^ 1 + u, ' 6 6 

+ sin^) ^ 1 t , ^-.-n + iA g ;e e 

(C.ll) 

The last four terms have e~ a 10 e _2 ^ T_T1 ^ as the argument of 3F2. This can be replaced simply 
by e~ 2< - r ~ ri - ) , and we obtain (13.431) . On the other hand, e~ CT1 ° for the first four terms cannot 
always be replaced by 1, since 3F2 of argument unity is singular if d — 3 — n > 0. 
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